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1. Motivations for Long-Range Modifications of Gravity 

Historically, the first attempt to modify the laws of gravity commonly accepted at that time was 
due to Laplace [jl]] who, in 1805, tried to add a velocity-dependent term to the standard inverse- 
square law of Newton to account for the finite velocity of propagation of gravity. But this work 
did not find echo practically until the surroundings of 1 880, when a series of works to estimate the 
gravitational finite propagation speed began. Such attempts to find deviations from the Newtonian 
inverse-square law of gravitation were performed to explain the anomalous secular precession of 
Mercury's perihelion, discovered by Le Verrier [^, without invoking undetected (baryonic) matter 
like the hypothesized planet Vulcan: for example. Hall [^] noted that he could account for Mer- 
cury's precession if the law of gravity, instead of falling off as 1/r^, actually falls of as \/v^ with 
k = 2.00000016. However, such an idea was not found to be very appealing, since it conflicts with 
basic conservation laws, e.g.. Gauss Law, unless one also postulates a correspondingly modified 
metric for space. Other historical attempts to modify Newton's law of gravitation to account for the 
Mercury's perihelion behavior yielded velocity-dependent additional terms: for a review of them 
see Ref. [Q] and references therein. Such attempts practically ceased after the successful explana- 
tion of the perihelion rate of Mercury by Einstein [||] in terms of his tensorial general theory of 
relativity: an exception is represented by Manev [^] who, with a 1/r^ correction to the Newtonian 
potential, was able to reproduce the anomalous apsidal precession of Mercury. 

Moving to more recent times, in the modern framework of the challenge of unifying gravity 
with the other three fundamental interactions of Nature, it was realized that possible new phenom- 
ena could show up just as deviations from the Newtonian inverse-square law of gravitation [^, |8|]. 
In general, they would occur at submillimeter length scales. Concerning general relativity, tradi- 
tionally, corrections to it, in the form of modifications of the Einstein-Hilbert action in order to 
include higher-order curvature invariants with respect to the Ricci scalar, were considered to be 
important only at scales close to the Planck one and, consequently, in the early universe or near 
black hole singularities^ ; see, e.g., Ref and references therein. It was not expected that such cor- 
rections could affect the gravitational phenomenology at low energies, and consequently at larger, 
macroscopic scales. 

Thus, why considering modifications of the standard laws of gravity at large, astronomical, 
astrophysical and cosmological as well, as done in recent years? To try to accommodate some 
recently observed phenomena, occurring at very different scales ranging from solar system to cos- 
mological distances, which, at present, have not yet found fully satisfactorily mundane explanations 
in terms of conventional physics, gravitational or not [ p^ [TT| , p^ . Below we list just some of such 
anomalous effects. 



• The flyby anomaly. It consists of unexplained changes of the asymptotic outgoing velocities 
of some spacecraft (Galileo, NEAR, Cassini, and MESSENGER) that occurred at their clos- 



est approaches with the Earth [13, 14]. Is it due to conventional non-conservative effects, or 



are modifications of the laws of gravity responsible of it? 



'it was shown that in such ways the non-renormalizability of general relativity became more tractable in the context 
of quantum field theory. 
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• The anomalous perihelion precession of Saturn [15] detected by processing almost one cen- 
tury of planetary observations with the inclusion of the latest radiometric data of the Cassini 



spacecraft [ ]16[ ] : is it a data processing artifact or a genuine physical effect? 



10 ^^-2 



The Pioneer anomaly. It is an unexplained acceleration Apio = (8.74 it 1.33) x 10 
approximately directed towards the Sun affecting the Pioneer 10/11 spacecraft after they 



passed the 20 AU threshold in the solar system [ ]17| , |18[ ]. Is it induced by some mundane 
non-conservative effects, or is it a sign of modifications of the laws of gravity? 

The dark matter issue. In many astrophysical systems like, e.g., spiral galaxies and clusters 
of galaxies a discrepancy between the observed kinematics of their exterior parts and the 
predicted one on the basis of the Newtonian dynamics and the matter detected from the emit- 
ted electromagnetic radiation (visible stars and gas clouds) was present since the pioneering 
studies by Zwicky [ |l9| ] on the Coma cluster, and by Kahn and Woltjer [20], Bosma [21], and 
Rubin and coworkers [22, U] on individual galaxies. More precisely, such an effect shows 
up in the galactic velocity rotation curves [|^] whose typical pattern after a few kpc from the 
center differs from the Keplerian 1 / y/r fall-off expected from the usual dynamics applied to 



the electromagnetically-observed matter. Does the cause of such a phenomenology reside in 
the action of still undetected (non-baryonic) dark matter whose dynamics is governed by the 
usual Newtonian laws of gravitation, or have they to be modified? 

The dark energy issue. In recent years, an increasing amount of observational evidence has 
accumulated pointing toward the fact that the universe has entered a phase of accelerating ex- 
pansion. Some of such observations are of direct, geometrical nature: standard candles like 



the supernovae Snia [ |25| , Pq ], gamma ray bursts [ [27| ] and standard rulers like the CMB sound 
horizon [^, |^]. Other ones are of dynamical nature like the rate of growth of cosmologi- 



cal perturbations [ ]30| ] probed by the redshift distortion factor or by weak lensing ||3 1|]. All 
these observations are converging towards a confirmation of the accelerating expansion of 
the universe, assumed homogeneous. They are successfully fitted by the simplest cosmolog- 
ical model predicting accelerating cosmic expansion: its ingredients are the the assumptions 
of flatness, validity of general relativity, the presence of a cosmological constant A, identi- 
fied with an unknown and still directly undetected form of energy (named dark energy for 
these reasons), and Cold Dark Matter (ACDM). However, for some puzzles of the ACDM 



cosmology, see Ref. [32]. Contrary to the ACDM paradigm, is it possible to accommodate 
the aforementioned observations by invoking modifications of the standard laws of gravity? 

Some of the models of modified gravity that have been proposed to address the aforementioned 
phenomenology are listed below. 



Dvali-Gabadadze-Porrati (DGP) braneworld model [|33|]. In it our universe is a (3-(-l) space- 
time brane embedded in a five-dimensional Minkowskian bulk. All the particles and fields 
of our experience are constrained to remain on the brane apart from gravity which is free 
to explore the empty bulk. Beyond a certain threshold ro, which is a free-parameter of the 
theory and is fixed by observations to about 5 Gpc, gravity experiences strong modifications 
with respect to the usual four-dimensional Newton-Einstein picture: they allow to explain 
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the observed acceleration of the expansion of the Universe without resorting to the concept 
of dark energy. 



Modified Gravity (MOG) by Moffat |34]. It is a fully covariant theory of gravity which is 
based on the existence of a massive vector field coupled universally to matter. MOG yields 
a Yukawa-like modification of gravity with three constants which, in the most general case, 
are running; they are present in the theory's action as scalar fields which represent the gravi- 
tational constant, the vector field coupling constant, and the vector field mass. It has used to 
successfully describe various observational phenomena on astrophysical and cosmological 



scales without resorting to dark matter [|35|]. 



Modified Newtonian Dynamcs (MOND) [^, ^ |38|]. It is a non-linear theory of gravity 
which predicts departures from the standai^d Newtonian inverse-square law at a characteristic 
acceleration scale [ p9[ ] Aq = 1.27 x 10^^" m s^^ below which the gravitational acceleration 
gets a « 1 /r behavior. MOND was proposed to explain certain features of the motion of 
ordinary electromagnetically detectable matter in galaxies and of galaxies in galactic clusters 
without resorting to exotic forms of still undetected dark matter. 



f{R) models [40]. These theories come about by a straightforward generalization of the 
Lagrangian in the Einstein-Hilbert action in which the Ricci scalar R is replaced by a general 
analytical function f{R) of R. They have mainly been used in cosmological and astrophysics 



scenarios [41] 



• Curvature Invariants models. They encompass inverse powers [^, 43] and logarithm [|4| of 
some curvature invariants in the Einstein-Hilbert action and have been used for tackling the 
dark energy-dark matter problems. 

• Yukawa-like models. There are many theoretical frameworks yielding such a modification 
of the Newtonian inverse-square law [45, 12]. Models encompassing Yukawa-type extra- 
accelerations have been used for a variety of applications ranging from solar system effects 



like the Pioneer anomaly [34] to astrophysical and cosmological scenarios [47, 48, 49, 5C, 



51, 52]. 



Hooke-like models. With such a definition we mean models of gravity introducing an addi- 
tional term proportional to the distance r. An important case is given by the Schwarzschild-de 



Sitter spacetime [ 53 ] which yields a correction to the Newtonian inverse-square law propor- 
tional to Ar [p4|], where A is an uniform cosmological constant. Another example is given by 
the extra-acceleration proposed by Jaekel and Reynaud [ ^5| ] to explain the Pioneer anomaly. 



Pioneer-like models. The simplest one consists of postulating a constant and uniform accel- 
eration radially directed towards the Sun and having the same magnitude of Apio existing in 
the outer regions of the solar system at heliocentric distances r > 20 AU. Other forms have 

5 



been postulated for it, both distance-dependent [pa, 56] and velocity-dependent [57 



■ 



Other models of modified gravity that we will not consider here aie the Einstein-Aether theory 
[]6T1], Tensor- Vector-Scalar (TeVeS) [52] and braneworld gravity [^]. 
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Table 1: Corrections Afjj, in 10^^ arcsec cy^' , to the standard Newton/Einstein perihelion precessions of 
the inner planets estimated by E.V. Pitjeva with the EPM2004 [|^ (Mercury, Eaith, Mars) and EPM2006 
(Venus) ephemerides. The result for Venus has been obtained by recently processing radiometric data 
from Magellan spacecraft (E.V. Pitjeva, personal communication, 2008). The errors in square brackets are 
the formal ones: the other ones have been re-scaled by Pitjeva to get realistic estimates. 



Mercury 


Venus 


Earth 


Mars 


-36 ±50 [42] 


-4±5[1] 


-2±4 [1] 


1±5[1] 



The only motivations for the aforementioned models are just the phenomena themselves for 
which they have been introduced. Thus, such models must be put on the test in different scenarios 
by devising independent observational checks. In particular, they must not exhibit discrepancies 
with the well-tested standard laws of gravity in local astronomical systems like, e.g., our solar 
system and the double pulsar PSR J0737-3039 [||, H]. 

2. Our Method for Testing Modified Models of Gravity in Astronomical Scenarios 

In general, a given Long Range Modified Model of Gravity (LRMMG) yields predictions 
for certain observable effects of the form 

^ = Kg, (2.1) 

where ^ ^ implies no modifications of gravity, and g is a function of the geometrical configu- 
ration of the system adopted characteristic of the LRMMG considered; the fact that the LRMMG 
parameter K enters as a multiplicative factor in the predicted effects ^ will be very important for 
us, as we will see below. For example, for the solar system's planetary longitudes of the perihelia 
CJ, it turns out 

g = g{a,e) (2.2) 

with a semimajor axis and e eccentricity of the planetary orbit considered. Corrections ACT to the 
usual Newtonian/Einsteinian perihelion precessions have been recently estimated for several plan- 



ets of the solar system by E.V. Pitjeva [pq, 67, Iq] by fitting 100 yr of observations of several kinds 



with the force models of various versions of the EPM ephemerides (EPM2004 Q] EPM2006 [g^], 
EPM2008 [p^]). Since they do not include any LRMMGs, the corrections AdT, by construction, 
account for them, so that they will be our : they are listed in Table |l]. By directly comparing AdJ 
to the predicted anomalous perihelion precession = Kg{a,e) for each planet separately allows 
to put upper bounds on \K\ since ACT are compatible with zero, according to Table |l|. This approach 
is good for test the hypothesis K = 0. The hypothesis K y^O can be tested by taking the ratios 

Uab = (2.3) 

of AdJ for different pairs of planets A and B, and comparing them to the predicted ratios 

_ g{aA,eA) . 
SAB = = —f r ■ (2.4) 
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Table 2: Ratios Hab of the estimated corrections of the perihelia AU5 for the pairs of planets A B. The 
uncertainties SHab have been evaluated by using the realistic errors for the individual corrections of Table 
I- 



A 


B 






Mercury 


Venus 


9 


24 


Venus 


Mercury 


0.1 


0.3 


Mercury 


Earth 


18 


61 


Earth 


Mercury 


0.06 


0.19 


Mercury 


Mars 


-36 


230 


Mars 


Mercury 


-0.03 


0.18 


Venus 


Earth 


2 


6.5 


Earth 


Venus 


0.5 


1.6 


Venus 


Mars 


-4 


25 


Mars 


Venus 


-0.2 


1.6 


Earth 


Mars 


-2 


14 


Mars 


Earth 


-0.5 


3.5 



i^AB> by construction, does not expUcitly depend on K, but it still retains a functional dependence 
on a and e typical of the LRMMG considered. If 

Hab / ^AB (2.5) 

within the observational errors, i.e. if 

»FAB = |nAB-^AB|/0, (2.6) 

the LRMMG considered is severely challenged. Of course, the uncertainty in the ratios ITab has to 
be evaluated in a realistic and conservative way to reduce the risk of unsound conclusions; see the 
discussion in Section Q We quote the ratios IIab for all the pairs of the inner planets in Table ^. 

An analogous approach can be followed with the double pulsar PSR J0737-3039. Now, one 
observable is Aft), i.e. the difference between the phenomenologically determined periastron pre- 
cession ftJmeas ^ud the Computed rate ftJipN for the general relativistic IPN rate (it is analogous to the 
well-known Mercury's perihelion precession by Einstein). Another observable is AP, i.e. the dif- 
ference between the phenomenologically determined orbital period Pb and the computed Keplerian 
one P^'^P. Thus, the observation-based ratio 

(2.7) 

AP 

can be compared with the predictions ^ for the same ratio by various LRMMGs; again, if they are 
equal within the errors, LRRMMG passes the test, otherwise it is challenged. 
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3. Results from Local Systems 

3.1 Dvali-Gabadadze-Porrati (DGP) Braneworld Model 



The DGP braneworld model []33[], put forth to explain the cosmic acceleration without resort- 
ing to dark energy, predicts an anomalous perihelion precession independent of the semimajor axis 

ft^DGP = =F — = =f0.0005 arcsec cy"\ ro 5 Gpc, (3.1) 
8ro 

where the =p sign are related to the two different cosmological branches of the model. Thus, 

^AB = 1 (3.2) 
for it. By linearly adding the absolute values of the uncertainties in AttJ, it turns out for the pairs 



A=Mars, B=Mercury and A=Earth, B=Mercury [72| 



^I'Mars Mercury = 1 .0 it 0.2, ruled out at 5a Icvcl (3.3) 

*PEarth Mercury = 0.9 ± 0.2, rulcd out at 4.5a Icvcl (3.4) 
In the case of the double pulsar [|73|], the predicted ratio ^dgp between CJdgp and /bo? is 

^OQP = (1.4±0.1) X lO"'' s-2. (3.5) 

Instead, the observed ratio Moh& between Aw = OJmeas — Wipn and AP = — P^'^^ is 

^n,eas = (0.3±4)xl0-'^ S-2. (3.6) 

They are not compatible at 14a level. 
3.2 Modified Gravity (MOG) 



MOG 1 34] predicts a Yukawa-type extra-acceleration [74, 75] 



GM , , . . 
Amog = ^a[l-(l+Atr)exp(-/ir)], (3.7) 

where a and pL have been fitted to a set of astrophysical galactic data [^] in the framework of the 
searches for an explanation of the flat rotation curves of galaxies without resorting to dark matter. 
For a Sun-planet system we have a perihelion extra-rate ( 1 //i = 33, 000 AU) [ [7^ ] 



G^mog ~ -WM YGMa(l -^2). (3.8) 

The inner planets yield [ f76| ] 

^I'venus Mercury = 1.3 ±0.3, rulcd OUt at 4a IcVCl (3.9) 

^J'Earth Mercury = 1.6 it 0.2, rulcd OUt at 8a IcVCl (3.10) 

^J'Mars Mercury = 2 ±0.2, rulcd OUt at 10a IcVCl (3.11) 
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3.3 Modified Newtonian Dynamics (MOND) 



MOND ^ ^ predicts acceleration-dependent modifications of the Newtonian inverse- 
square law to explain the almost flat galactic rotation curves 



. ^Newton ^ 



(3.12) 



where galaxies data fitting yields [ 39] Aq = 1 .2 x 10 m s ^ and (x) is an interpolating function 
whose most widely used forms are 



|77], M 



l+x 



|78 



with 



and 



;U ^ 1, x> 1, 



H ^ X, X <^ \. 



(3.13) 



(3.14) 



(3.15) 



solar system observations are, strictly speaking, tests of the form of the interpolating function in 
the large acceleration limit (x » 1). 

The anomalous perihelion precession predicted by MOND is, in the large acceleration limit 



(x> 1), [36, 7 



C^OND = —kon 
^0 = 1, h = I correspond to 

ko = 1/2, h = 2 correspond to 



''MOND 



2ft 



h, n 



GM 



- , ''MOND 



GM 



l+x' 



(3.16) 



(3.17) 



(3.18) 



It turns out [ pT] ] that the ratios ITab of the perihelia, independent by construction of ko, for the 
pairs A=Mars, B=Mercury and A=Earth, B=Mercury rule out the MOND perihelion precessions 
for 1 < /j < 2 at several sigma level. 

3.4 Extended /(/?) theories 



Extended theories of gravity with [ ]82p 
where R is the Ricci scalar, yield power-law corrections to the Newton's law 



(/3-l)GM 



(3.19) 



(3.20) 
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which, for j3 = 0.817, obtained by fitting the SNela Hubble diagram, yields good results in fitting 
some galactic rotation curves; j8 is related to the exponent k of the Ricci scalar R. The induced 



perihelion precession is [ ]83| ] 



0^m= ^^~^, -^-"'G{e,^)^ (3-21) 

where G{e,p) is a function of the eccentricity and j3. The resulting ^ab are not compatible with 
zero for many pairs A and B of inner and outer planets at several sigma level [^3|]. 

3.5 Inverse powers and logarithm of some curvature invariants models 

Inverse powers of curvature invariants in the action lead to [^, 43] 



aGMi6k + 3) 



For k = I and the Sun, fc = 10 pc, so that 



45 a 

m=i ~ ^\/GMa'T{l-e^) (3.23) 



' c 



Thus, the inner planets yield [ [7^ ] 

^I'Mars Mercury = 10^ zb 0. 1 , ruled OUt at 10^0 IcVCl (3.24) 

^Mam Earth = 38 it 3.5, rulcd out at 11a level (3.25) 

*I'Earth Mercury = 10^ ± 0.2, ruled out at lO'^a Icvcl (3.26) 
Logarithm of some curvature invariants in the action, able to modify gravity at MOND scales 



so to jointly treat Dark Matter and Dark Energy, induces [44] 

^2 



ALog- r, = 0.04pc, (3.27) 



so that 



G^Log - ^ (3.28) 



The pairs A=Mars, B=Mercury and A=Earth, B=Mercury yield [72] 



*J'Mars Mercury = 30.7ib0.1, rulcd out at 31a Icvcl (3.29) 
^I'Earth Mercury = 10.6 it 0.2, rulcd out at 53a Icvcl (3.30) 
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3.6 Yukawa-like theories 

Many theoretical frameworks ^,12] yield a Yukawa-type acceleration 



GMa / , r\ f r . ^ . ^ 

AYukawa = ^ + X/ '^'^ V X ' ' 



For A > ae, the induced perihelion rate is [ 84 ] 



aVGMa ( a 
C^Yuk ~ —ZTi — exp - - , (3.32) 



2A2 ^'^^V A 
so that for a pair of planet A and B 

The pair A=Earth, B=Mercury allows to obtain [ p4| ] 

A =0.182 ±0.183 AU. (3.34) 



Solving for a and using AdJ and a for Mars yields [ 84 ] 



a = ?^^exp ( ^) = (0.2 ± 1) x 10-^ (3.35) 



in which we have used the value for A obtained from Earth and Mercury. 

3.7 Hooke-type theories 

A Hooke-like acceleration is induced, e.g., by a uniform cosmological constant A in the 



Schwarzschild-De Sitter spacetime [54] 



Aa = ^AcV. (3.36) 



The induced perihelion shift is [85] 



Some pairs of planets yield 

^I'Mars Mercury = 7.8 ± 0.2, rulcd OUt at 39(7 IcVCl, (3.38) 
^I'Earth Mercury =4.1 ±0.2, rulcd OUt at 20a IcVCl, (3.39) 
^Jupiter Mercury = 51 ± 12, rulcd OUt at 4a IcVCl (3.40) 

In the case of the double pulsar, the predicted ratio between Wa and Pa is [f73[| 



^A = (3.4 ± 0.3) X 10"^ s-^. (3.41) 
Instead, the observed ratio ^ohs between Ad) = ci)meas — ^ipn and AP = P^ — P^'^ is [73] 

i^n,eas = (0.3±4)xl0-" S-2. (3.42) 

They are not compatible at 1 1 a level. 
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3.8 The Pioneer anomaly 

If the Pioneer Anomaly ||17, 18] was of gravitational origin, the exotic force causing it should 
also act on the planets of the Solar System, at least on those moving in the spatial regions in which 
it manifested itself in its presently known form. A constant and uniform radial acceleration with the 



same magnitude of that causing the Pioneer Anomaly would induce orbital effects too large 087| ] to 



have escaped so far detection, as shown by the AdT estimated with the EPM2006 ephemerides [ |670 



for Jupiter, Saturn and Uranus []88|,|89|], and from the residuals of right ascension a and declination 
8 computed with the JPL ephemerides [|90|]. Also some recently proposed velocity-dependent 
exotic forces 1 57, 58, 6^ are ruled out by planetary observations [90,^]. For other works dealing 



with the same problem, see Refs. [92, 93, 94, 95]. 



3.9 The anomalous perihelion precession of Saturn 

E.V. Pitjeva has recently fitted almost one century of planetary data of all kinds, including also 



3 years of radiotechnical data from Cassini, with the EPM2008 ephemerides [[Iq]. She estimated 
(E.V. Pitjeva, private communication, 2008) a statistically significant non-zero correction 



-0.006 ±0.002 arcsec cy" 



(3.43) 



while the formal error of the fit is 0.0007 arcsec cy^' . However, applying the standard re-scaling 
by a factor 10 of the optical only observations there is also, in principle, the possibility that the 
uncertainty can be as large as 0.007 arcsec cy^^ (E.V. Pitjeva, private communication, 2008). Pre- 
vious results obtained with the EPM2006 ephemerides ||67|], which did not include Cassini data, 
yielded 

AGTsatum = -0.92 ± 0.29 arcsec cy"^ (formal error). (3.44) 
Is it really a genuine physical effect needing explanation, or is it some artifact of the data reduction 



procedure? None of the exotic models examined so far is able to accommodate it [15|. 



3.10 The general relativistic Lense-Thirring effect 

Until now we have only dealt with putative modifications of the standard Newtonian/Einsteinian 
laws of gravity. In fact, the estimated corrections AdJ to the usual rates of the perihelia account, 
by construction, also for a standard general relativistic effect which has not been included in the 
dynamical force models of the EPM ephemerides used to determine them: it is the gravitomagentic 
Lense-Thirring effect [ p6| ] consisting of secular precessions of the form 

2G5(l-3cos/) 



c2a3^1-^'2)3/2 ' 



(3.45) 



where S is the Sun's angular momentum and / is the inclination of the planetary orbital plane to the 
Sun's equator which is quite small for all the inner planets. The Lense-Thirring precessions for the 
inner planets are of the order of 10"^ _ iq-^ 

arcsec cy ^ , and their direct measurability has been 



discussed in, e.g., Ref. []97[]. It is interesting to check if the unmodelled Lense-Thirring precessions 
pass the test of the ratio of the perihelia, with 



Sab 



4(1 



4(1 



-.2)3/2- 



(3.46) 
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Table 3: Ratios of the predicted Lense-Thirring precessions for the pairs of planets A B. 



A 


B 


i:LT 
bAB 


Mercury 


Venus 


7 


Venus 


Mercury 


0.1 


Mercury 


Earth 


18 


Earth 


Mercury 


0.05 


Mercury 


Mars 


64 


Mars 


Mercury 


0.01 


Venus 


Earth 


3 


Earth 


Venus 


0.4 


Venus 


Mars 


9 


Mars 


Venus 


0.1 


Earth 


Mars 


3.5 


Mars 


Earth 


0.3 



By comparing Table ^ with Table ^, it can be seen that *Pab is compatible with zero for all the pairs 
A B of inner planets, contrary to all the models of modified gravity examined so far. 

4. Discussion and conclusions 

4.1 Some technical aspects 

If and when other teams of astronomers will independently estimate their own corrections ACJ 
with different ephemerides, it will be possible to fruitfully repeat all the tests presented here. In 
doing them we always used the formal errors in Aijj re-scaled by a factor 2 — 5 for the inner 
planets and up to 10 times for the outer planets for which mainly optical data have been used. 
Moreover, in view of the correlations among the estimated Adl, we always linearly propagated 
their errors, instead of summing them in quadrature, by getting 

SHab < Hab + . (4.1) 

For example, the correlations between the corrections AdJ of Mercury and the Earth are as large 
as 20% (Pitjeva, private communication, 2005). Concerning the uncertainties in <^ab, in principle, 
they should have been computed by propagating the errors in the semimajor axes a and the eccen- 
tricities e of the planets A and B entering them for each model considered. However, they are quite 
negligible because the relative (formal) uncertainties in the semimajor axes and eccentricities of 
the inner planets are all of the order of 10^^^ and 10^^ — 10^' \ respectively, according to Table 3 
of Ref. []67[]. Processing more ranging data from Cassini will help in clarifying if the perihelion 
precession of Saturn can really be considered as a genuine physical effect. 
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4.2 Conclusions 

• All the long-range modified models of gravity examined here are severely challenged either 
by the ratios of the periheha of different pairs of solar system's planets or by the double 
pulsar combined data for the periastron and the orbital period. Both such kinds of ratios 
cancel out the small multiplicative parameters which directly account for the various exotic 
models considered, but they still retain a pattern which is characteristic of the model tested 
through the orbital parameters of the systems used. Only the general relativistic Lense- 
Thirring effect, not modelled in the EPM ephemerides and, thus, accounted for, in principle, 
by the estimated corrections ACJ, passes the test of the ratio of the perihelia. 

• None of the exotic models examined is able to explain the anomalous perihelion precession 
of Saturn. The analysis of more Cassini data will help in clarifying if it is really incompatible 
with zero at some statistically significant level, thus requiring an explanation in terms of some 
physical phenomena, or if it is some artifact of the data reduction procedure. 
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